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Abstract 

We develop a conditional sampling scheme for pricing knock-out barrier options under the Linear 
Transformations (LT) algorithm from Imai and Tan (2006), ref. [12]. We compare our new method 
to an existing conditional Monte Carlo scheme from Glasserman and Staum (2001), ref. [9], and show 
that a substantial variance reduction is achieved. We extend the method to allow pricing knock-in 
barrier options and introduce a root-finding method to obtain a further variance reduction. The 
effectiveness of the new method is supported by numerical results. 

1 Introduction 

Barrier options are financial instruments whose payoff depend on an underlying asset hitting a specified 
level or not during some period. The simplest example would be a regular knock-out call option, where 
the payoff is set to zero if the option hits a certain barrier level during the lifetime of the option, and 
results in the payoff of a regular call option otherwise. These options are popular because they cost at 
most the amount of their vanilla counterparts as there is a higher probability of ending up with nothing. 
For more background information on barrier options and their applications, see for instance [2], [3], [24] 
and [25]. For sufficiently simple products there exist analytical formulas for the value, see, e.g., all the 
previous references, or [11]. When the product is too complex, one must resort to numerical methods to 
obtain the value. Furthermore, if the barrier is not monitored continuously but discretely, no accurate 
closed form formulas are available. Several papers exist dealing with discretely monitored barrier options, 
e.g., [14] and [26]. 

A straightforward simulation method for pricing discretely monitored options is to sample the assets 
in each monitoring date. However, when a barrier condition is introduced in the payoff, only a minority 
of the paths may lead to non-zero values if the probability of satisfying the barrier condition is low. 
This leads to significant increases in the variance of the estimator. Glasserman and Staum [9] remedy 
this problem by introducing a sampling scheme where the asset with the barrier condition is sampled 
conditional on survival in each monitoring date. In [14] a different scheme is used based on the hitting 
time of the asset as the sampling variable instead of the asset increments. Both of these papers use 
Monte Carlo to sample the necessary variates. 

In this paper we use low-discrepancy points to generate the uniform variates for sampling the asset 
paths. The merit of using low-discrepancy points for valuing financial products is discussed in [21] and 
[22]. Since then quasi-Monte Carlo methods have often been applied to applications in finance, e.g., 
[16, 1, 7, 15, 18]. Care must be taken however in the construction of the asset paths, see for instance [19] 
or [23] , where it is shown how the path construction method can have a large influence on the convergence 
rate. Imai and Tan [12] provide an alternative to the Brownian Bridge or PCA construction, called the 
Linear Transformation method (or LT for short). They provide an optimal construction, in some sense, 
by minimizing the effective dimension of the underlying payoff function. We reconcile this LT method 
with the idea of conditional sampling, to gain a "best of two worlds" approach for valuing barrier options. 

The outline of our paper is as follows. In Sections 2 and 3 we present an overview of the LT algorithm 
from [12] and the conditional sampling scheme from [9], Section 4 contains the main results of our paper: 
the conditional sampling scheme for the LT algorithm. In Section 5 we give several numerical examples 
to illustrate our method. Then in Sections 6 and 7 we present two extensions to our sampling scheme. 
The first extension deals with knock-in options, and the second one is a modification to combine the 



quasi-Monte Carlo method with an analytical integration using root-finding. We summarize our results 
in Section 8. 



2 The LT method 



We assume a Black-Scholes world in which the risk-neutral dynamics of the assets are given by 

dSi(t) = rSi(t)dt + <TiSi(t)dWi(t), i = l,...,n, 

where Si(t) denotes the price of asset i at time t, r is the risk- free interest rate and <Ji the volatility of 
asset i. Also, XV = (Wi(t), . . . ,W n (t)) is an n-dimensional Brownian motion, with dWidWj — Pijdt. 
For more information on this type of market model, see for instance [5], [8], [11] or [24]. When resorting 
to Monte Carlo techniques for pricing options under this model, asset paths need to be sampled at each 
time point. In what follows, we assume an equally spaced time discretization At — T/m for simplicity, 
but we note that all results can be applied without using this assumption. Under this assumption, we 
have tj = jAt. Under Black-Scholes dynamics we simulate the trajectories according to 

s t {t) = s l {Qy {r ^ 2j2)t+a ' m{t) - 

Write XV = (oiWi(ii), . . . , GiW\{t m ), c^Wa^i), ■ ■ ■ , cr n W n (t m ))' , where we use the prime to denote the 
transpose of a vector. Then XV is multivariate normally distributed with covariance matrix 
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One approach to simulate XV is to decompose the matrix £ as 

CC" = S 



where C is the Cholesky factor (see, e.g., [8]) and calculate 

w = Cz (1) 

where z is an mn-dimensional vector of i.i.d. standard normally distributed variables. Assuming that 
we have a European option payoff represented as 

max 0) 

the standard Monte Carlo method simulates the function f(W) by mapping a uniform variate u in the 
unit cube to z by applying the inverse cumulative distribution function When using quasi-Monte 

Carlo, the uniform variates are replaced by a low-discrepancy point set. The matrix C in (1) can be 
multiplied with any orthogonal matrix Q while retaining the correct distribution for w, see, e.g., [8]. 
Under the LT method described in [12], instead of simulating / directly from z, first an orthogonal 
transformation is applied, i.e., / is simulated from 

W = Az = CQz (2) 
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for a carefully chosen orthogonal matrix Q. In what follows we use the notation Q.k to denote the 
fcth column of Q and Qk- to denote the fcth row. The matrix Q is chosen according to the following 
optimization problem: 

maximize variance contribution of / due to kth dimension 

subject to IIQ-fcll = 1, 

(Q*,Q. fe )=o, j = i,...,*-i, 

where Q*- denotes the columns of Q that have already been optimized in the previous iterations. The 
algorithm is carried out iteratively for k = 1, 2, . . . , mn so that in the fcth optimization step the objective 
function ensures that, given columns Q*-, j — 1, . . . , h— 1 which have already been optimally determined 
in the previous iterations, the variance contribution due to the fcth dimension is maximized while the 
constraints ensure orthogonality. Being able to express the variance contribution for each component 
analytically for general payoff functions / can be quite complicated. Therefore, Imai and Tan [12] propose 
to approximate the objective function by linearizing it using a first-order Taylor expansion around a point 
z = z + Az, 

Az e . 

Using this expansion, the variance contributed to the kth component is 




The expansion points are chosen as Zf. = (1, . . . , 1, 0, . . . , 0), the vector with k — 1 leading ones. The 
optimization problem becomes 

maximize 
subject to HQ-fcll = 1, 

(Q*pQ- k ) = o, .i i k 1- 

In [12], only the first 25 or 50 columns are computed. The loss of efficiency is argued to be minimal, 
particularly when the underlying LT construction is effective at dimension reduction with the first few 
dimensions already capturing most of the variance. However, as this construction is just part of the 
startup cost, we construct the entire matrix in our numerical tests. 



libit f~\ £ 




The original paper, [12], considers a basket arithmetic average option (Asian basket) to illustrate the 
computational advantage of the LT method. The payoff function inside the max-function in this case is 

mn 

f(z) = e Mi+E ^" a ^ - K 

i=l 

where K is the strike price, and 



[M = \og(w lui2 S ll (0)) + \ r - t i2 , 

where j 2 corresponds to the weight given to asset ii ~ — \)/m\ + 1 at time %2 = i — («i — l)m. 
From the first-order Taylor expansion we have 

mn / mn \ 



e=i \i=i 
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The optimization problem becomes 



\ 2 

fc— 1/ 



maximize V e ^+TC<^-^> (C,, Q. k ) 

subject to IQ-fel = 1, 

(Qr,-,Q.fc)=Q, i = i,...,fc-i. 



For fc = 1, the solution is found as 
for subsequent fc the solution is 



gie ME J t 1 1 (a^) (C . )Q , fc) 
|| E - e M i+ E^(c, ! ^> (c . )Q ^ )|| ' 



where the sign can be chosen freely. The proof is in the appendix of [12]. Later in Section 4 we revisit 
the sign choice. 



3 Conditional sampling 

Suppose we are interested in pricing an up-&-out option on n assets with payoff 

0(Si(ti), ■ • • , S n (t m )) = max (/(5i(ti), . . . , S n (t m )), 0) I jmax Si(tj) < B 

where I is the indicator function. The valuation can be done by sampling the assets using Monte Carlo 
and then calculating the payoff g(S\(t\), . . . , S n (t m )). When the probability of survival becomes low 
(e.g., if B — Si(to) is small) a lot of the generated paths will result in a knock-out. This can make 
the variance among all paths quite large relative to the price. To remedy this problem Glasserman and 
Staum [9] propose an estimator using conditional sampling. Given S\(tj) for j = 0, . . . , i — 1, define 

i-l 

L t = l[¥[S 1 (t J+1 )<B\S 1 (t 1 )}. (4) 

3=0 

The variable Li can be interpreted as the "likelihood" of the path surviving i steps. Under the Black- 
Scholes framework, the probability in the product can be analytically determined as 

io ste)H r -^) At ~ 

The authors [9] now sample the variable 

,9i = L m max(/(5i(*i), . . . ,S n (t m )),0) , 

where all Si(tj) < B, j = 1, . . . , m. Therefore, there are no paths knocked out in the simulation, greatly 
reducing variance if such an event has a high probability. An overview of this algorithm is presented in 
Algorithm 1. We recapitulate two theorems from [9], which we prove using methods we will use again 
later. 

Theorem 1. The estimator based on sampling 

.9i = i m max(/(S , i(ii), . . . ,S„(i m )),0) (5) 
is unbiased. That is, E[gi] = E[g]. 



F[S 1 (t j+i )<B\S 1 (t j )] = $> 
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Proof. Denote by x the function expressing the barrier condition 



X(U) = 



1 if max ie{1 ... im} S(tj) < B, 
otherwise, 



where we explicitly denote the dependence of \ on U , the uniform variates used to sample the asset 
paths. To stress the dependency of / and Tj on U we will write this explicitly as well in terms of U . 
Furthermore, without loss of generality, we only consider the case of one asset and set r = 0. 



We find 



E [ X {U) max(/([/), 0)] = / X (u) max(/(«i, . . . , u m ), 0) du 

J[oa\™ 

= •••/ max(/(«i,...,u m ),0) du m ---dui. 

Jo Jo 

By using the change of variables u\/Y\ — ii\ to u m /T m (ii\, . . . , u m -i) = u m we obtain 



„ m 

E[ X (t/)max(/(J7),0)]= / T x TT T^, . . . ,r^( 



-l(ttl, . • .,Uj-l)) 

x max(/(r 1 u 1 , . . . , r m (ui, . . . , u m -i)u m ), 0) du\ ■ ■ ■ du Tl 



E 



L m max(/(C/),0)j , 

which is what we needed to prove. □ 

Furthermore, it can be proven that the variance of the new estimator is at most that of the original 
one. 

Theorem 2. The estimator using (5) has reduced variance. That is, Var[<7i] < Var[<?]. The inequality 
is strict i/ P [max-,- Si(tj) < B] > and E[/] > 0, i.e., if there is any chance of knock-out and positive 

payoff. 

Proof. Because of Theorem 1 we only need to prove that E[g1] < E[g 2 ]. We assume again that r = 
and n = 1. Using the same notation as in the proof of Theorem 1 we can write 



mi] 



/ l? m max(/(u), 0) 2 du 
J[0,ir 

X 2 (u) L m max(/(w), 0) 2 du 

[o,i]- 

< / x 2 (M)max(/(M),0) 2 du 
i[o,i]- 

where the inequality follows because < L m < 1 and the last equality follows since x 2 ( u ) = x( u )- D 



4 Conditional sampling under the LT method 

We now derive a conditional sampling scheme for the LT algorithm discussed in Section 2. As an example 
we will again consider the up-&-out option with payoff 

g(S 1 (t 1 ),...,S n (t m ))=max(f(S 1 (t 1 ),...,S n (t m )),0)llmaxSi(t j ) < B 
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Algorithm 1 The Glasserman and Staum algorithm [9] for an up-&-out option. 



Input: N, n, m, 5i(0), . . . ,S n (0), a x , . . . ,er n , (pij)i,j=i,..., n , r, qi, . . . ,q n , h, . . . ,t m = T, B, f 
Generate N mn-dimensional uniform vectors vM\ i = 1, . . . , N 
Calculate the matrix C 
Set L = 
for i = 1 to N do 
for j — 1 to m do 

^-Dn+l = P [Si(t J+1 ) < BlSifa)] v^_ 1)n+1 

Calculate Si(tj), . . . , S n (tj) 

Set L = L x¥[Si(t j+1 ) < B\Si(tj)} 
end for 

g l = Lx max{/(5i(ti), . . . , S n (t m )), 0} 
end for 

1 

° = \l nw^t) EjLi (ft - A) 2 



For the option to stay alive we have the restriction 



CTxWxitj) <log 
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5i(0) 



for all j = 1, 



.TO, 



(6) 



where under the LT method 

.7! Wife) = (A$" X H) . 

and the vector it consists of rrm uniform random variables. The restriction (6) becomes 



y2 a j,i® 1 ( u i) < kw f ° r a11 i = i. 



, to. 



We assume for the moment that a j i > for j = 1, . . . , m, which is true if Cj j > for j = 1, . . . , m 
and i = 1,. . . , uto and Q-,i appropriately chosen (cf. (2) and (3)), and later provide the solution for the 
general case. Using this assumption we can modify the inequalities such that there is only one constraint 
on the first uniform variable assuming the others are given: 



a i,i 



for all j = 1, 



Equivalently we can write these m conditions as the single condition 

"%) - aj j2 $ _1 ( u 2) - . 



U\ < $ I min 



a i,i 



=T„(ll2,...,«mn) 

Our algorithm samples u\ to M mn , then calculates the upper bound T u on u\ using u 2 , ■ ■ ■ , u mn , and 
then rescales u\ to Ui = Y u (u 2 , ■ ■ ■ ,u mn )ui in order to satisfy the barrier condition. Our estimator is 
thus based on sampling 

92 = T tl (M 2 ,...,u m „)max(/(Mi,u 2 ,...,u m „),0). 
For a down-&-out option, a similar analysis can be done, leading to the condition 

~b{tj) - a i:2 $ _1 (M 2 ) - ... - a^ mn $- l {u mn ) 



Li > $ ( 



max 



a j,i 



=Td{u 2 ,...,u mn ) 



G 



In case C has negative elements, it is possible to have negative values for a^\. Our algorithm can be 
modified to work around this problem: suppose a^i > for j 6 P C {1, . . . , m} and Oj.i < for all 
other j ^ V . Then we find the conditions on u\ as 

< b(t ] )-a^-Hu 2 )-...-a ] ^-\u mn ) for ^ g p 

and 

The restriction on now consists of an upper and a lower bound, denoted by T„ and respectively, 
suppressing the dependence on u 2 to it mn for ease of notation. The rescaled variable is then ui — 
Y c i + (T„ — Td)u\. The estimator is then based on sampling 

g 2 = max(T u - T d ,0) max(/(wi, m 2 , . • . ,u mn ),0). 

We take the maximum of T u — and as it could happen that T u < T d , and in this case the value of 
g 2 should be taken zero. Examples of such payoffs are given in Section 5.2 and Section 5.3. An overview 
of this algorithm is presented in Algorithm 2. To obtain a standard deviation we need to obtain a small 
number M of randomized estimators which is standard practice in applying QMC. Algorithm 2 is generic 
in the sense that one can either use a lattice rule or lattice sequence with random shifting (see [20] for a 
general reference) or a digital net or digital sequence with digital shifting (see [4] for a general reference). 
For our numerical tests we will use a digital sequence (being the Sobol' sequence with parameters from 
[13]) with digital shifting. 



Algorithm 2 The conditional LT algorithm 



Input: M, N, n, to, Si(0), . . . ,S n (0), cti,...,ct„, (pij)ij=i,..., n , r , ?i> • • • ,Qn, *i,---,*m = T,B,f 
Generate N TOn-dimensional low-discrepancy points u^', i = 1, ...,N 
Generate M random shifts AW, i = 1, . . . , M 
Construct the matrix A using the LT algorithm 
for k = 1 to M do 
for % — 1 to N do 

Randomize using and call this ■y( 2 ' fe ) 



Calculate the bounds and T u using v^'^ to »t'i? 
0i = T d + (T u - T d ) v[ l ' k) 

g kti = max(T u - T d , 0) max (/(A*- 1 ^.*)))^) 
end for 
end for 

~ _ l ^ w sr^M 
V ~ WW L^k=\ l~ii=\ 9k,i 



M(M 



2 

9k,i - A* J 



Similar as in the previous section we obtain the following results. 
Theorem 3. The estimator based on sampling 

g 2 = max(T u - T d ,0) max(/(«i,tt 2j • . . ,u mn ),0) (7) 
is unbiased. That is, E[<jf] = £[.92]- 
Proof. Denote by x the function 



1 if max ie {i ) ... im } S(ti) < B, 
otherwise. 
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Without loss of generality we assume r = 0. The value of the contract is then given as E [x{U) max(/(Z7), 0)] 
where we write out the dependencies explicitly in terms of U. We now want to prove that this equals 



E 



max(T„ - T d , 0) max(/(Z7), 0) 



E[x(l7)max(/(C7),0)] = / x( u ) max(/(tn, . . . , u mn ), 0) du 

J[0,1]™ 



max(/(ui, . . . ,u mn ),0) dui---du mn . 

[ ,l]m«-l J[T d ,T u ] 

By using the change of variables (ui — Td)/(T U — T^) = Ui and Ui = for i = 2, . . . , mn we find (using 
the notation u = (tti, U2, ■ ■ . , u mn )) that 

E[ X (C/)max(/(J7),0)] 

= / / max(T u - T d! 0)max(/(T (J + (T„ - T d )ui, . . . ,u mn ),0) dux ■ ■ ■ du mn 

J[0,l]mn-1 7 [ 0i l] 

= / max(T„ - Y d ,0)max(/(T d + (T u - T d )u 1: . . . ,u mn ),0) du 
= E [max(T„ - T d , 0) max(/(C7), 0) 

which is what we needed to prove. □ 

Theorem 4. When using regular Monte Carlo, the estimator defined in (8) has reduced variance. That 
is, Var < Var [g] . The inequality is strict ifV [max^ S±(tj) < B] > and E [/] > 0, i.e., if there is 
any chance of knock-out and positive payoff. 

Proof. Because of Theorem 3 we only need to prove that E[g|] < E[g|]. Using the same notation as in 
the proof of Theorem 3 we can write 

E\gl] = / max(T„ - T d , 0) 2 max(/(T rf + (T u - T d )ui, . . . ,u mn ),0) 2 du\ ■ ■ ■ du rnn 
V[o,i] m " 



[0,1]* 



X 2 (C/)max(T u - T d , 0) max(/(ui, . . . ,u m „),0) 2 dut ■■■du 7] 



< / x 2 (^ r )max(/(Mi,...,u mn ),0) 2 du x ■ ■ ■ du mn 

J[0,1]™ 

= E[ff 2 2 ], 

where the inequality follows because < < T u < 1 and the last equality since x 2 (U) — x(U)- □ 

For our case, when using quasi-Monte Carlo, we have the following theorem on variance. 

Theorem 5. When using a randomly shifted quasi-Monte Carlo method, the estimator defined in (8) 
has reduced variance. That is, Var[<72] < Var[<?]. The inequality is strict if P [max, , S± (tj ) < B] > and 
E [/] > 0, i.e., if there is any chance of knock-out and positive payoff. 

Proof. When using lattice rules or lattice sequences the standard randomization technique is to use 
random shifting [20] . For digital nets and digital sequences a standard technique is to use random digital 
shifting [4]. Results on variance for both these methods can be found in [16] and [17]. As random shifting 
is easier to understand we will use this in the proof, but the proof can be trivially modified for random 
digital shifting. In what follows we will use the notation {x} = x — \x\ to denote the fractional part of 
x. When using M > 1 random shifts the estimator becomes 



M N 

m E ^ E max ( T - - T <*' °) max (/ + ( T - - T >i*' fc) . v 2 ,k) . • • • . V ™) , o) 



k=l i=l 
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where the bounds T u and are calculated based upon the shifted points and these shifted points are 
obtained as 

v (%,k) _ | u (0 _|_ ^(fe) | for j = 1, ... , nm, 

with the random shift vectors A^ i.i.d. uniform variables over [0, l] nm . For ease of notation, write 

+ A«}) = / (T d + (T„ - T^U*'*', . .,^ n fc) ) • 

We will keep suppressing the arguments of the bounds and T u for ease of notation, but as a reminder 
we note that these depend on the integration variables as well. By observing that the integration of the 
linear map A i->- {u + A} over [0, l) mn is the same as integrating A itself over [0, l) m ™, we can easily 
prove that the above estimator is equivalent to E[§2]- Indeed, 



E 



M N 

M E jV E max ( T « - T <*< °) max(F({nW + A«}), 0) 

k=l i=l 
1 " 1 N 



E i E I max ( T « - T ^ °) max ( f ({« (,) + A(fc) »> o) dA(k) 

M N . 

E^E / max(T u ~T d ,0)max(^(A( fc )),0)dAW) 



M ^ N 

k=l i=l 
M N 

-Y- 

fe=i i=i- / l u - i ) 

M JV 

— V-VE[g 2 ] 
M N U2j 

fc=l i=l 



= E[ff 2 ]. 

Similarly, we find for the second moment 



E 



A 



1 

M 



M JV \ 

E E max ( T « - T ^ °) max(F({n« + A<"}), 0) 
fc=i i=i / 

M N \ " 

E ]v E max ( T « - T ^> °) max(F({u« + A«}), 0) dA« • • • dA< 

fe=i «=i / 

M JV \ 2 

E jv E max ( T « - T rf> °) max(F(A( fe )), 0) dA« ■ ■ ■ dA {M ^ 

fc=l i=l I 



I 


-7 ( 

J [0.1)™" \ 


/[0,l)m» 


-7 ( 

J[0,l) mn \ 



1 

M 



1 

M 



^ JVf-AT . 



M N 
M 2 N 

.u-.v 2 ' 

M^JV 



max(T n - T d , Of max(F(A), 0) 2 dA 



—2 / max(T u -T d ,0) X 2 (A)max(/(A),0) 2 dA 



which concludes the proof. □ 



5 Numerical Results 

In the numerical examples we use the Sobol' sequence with parameters from [13] and digital shifting 
[4] and we calculate all the columns of Q for the LT construction. We use this same setup also for the 
examples of the two extensions in Sections 6 and 7. 
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HAAr 1 I T T 1 I PQ 




(P, (Tl,P, A) 


OA /TO i T T l PC 


OA/TO I T T 


(Pi, 0. 25, 10000, 70) 


2682% 


26827o 


/ n r\ r»£ 1 AAAA TA\ 

(P2, 0.25, 10000, 70) 


257570 


25757o 


(Pi, 0.25, 125, 70) 


786% 


154% 


(P 2 ,0.25, 125,70) 


621% 


197% 


(Pi, 0.25, 105, 70) 


368% 


128% 


(P 2 , 0.25, 105,70) 


414% 


170% 


(Pi, 0.25, 110,100) 


231% 


120% 


(P 2 , 0.25, 110, 100) 


304% 


165% 


(Pi, 0.55, 105, 70) 


287% 


110% 


(P 2 , 0.55, 105,70) 


234% 


144% 


(Pi, 0.55, 105, 90) 


190% 


88% 


(P 2 , 0.55, 105,90) 


157% 


105% 


(Pi, 0.55, 150, 110) 


278% 


135% 


(P 2 , 0.55, 150, 110) 


269% 


126% 



Table 1: Single barrier Asian basket. The reported numbers are the standard deviations of the MC+CS 
method divided by those of the QMC+LT and QMC+LT+CS methods. The MC+CS method uses 
163840 samples, while the QMC+LT and QMC+LT+CS methods use 4096 samples and 40 independent 
shifts. 



5.1 Single barrier Asian basket 

Consider an Asian barrier option on four assets: 

(^ 4 130 

i=l j=l 

We will consider the valuation of this option under several model parameters. The fixed parameters are 
Si(0) — 100 for i = 1, . . . , 4, er 2 = 03 = 25%, 04 = 35%, r = 5% and T = 6 months. Here we have taken 
m = 130. We consider two correlation matrices: 

0.5 0.6 0.2 \ 
1 -0.2 -0.1 
0.2 1 0.25 ' 
0.1 0.25 1 / 

In what follows, we denote the LT method using a low-discrepancy point set by QMC+LT, and the 
conditional sampling scheme we derived for the LT method using such a point set by QMC+LT+CS. The 
conditional sampling scheme for Monte Carlo will be denoted by MC+CS. We report the improvement 
of the standard deviations of the value estimates of these methods to that of the conditional sampling 
scheme for regular Monte Carlo (MC+CS) for various choices of a\, K and B in Table 1. 

The first example takes o\ = 25%, K = 70 and the extreme choice of B = 10000. Clearly, the barrier 
will almost never be hit, and we should obtain the same result for QMC+LT as for QMC+LT+CS. This 
example also shows the power of the QMC+LT method compared to using MC+CS, as the standard 
deviation under the former method is about 26 times smaller. 

We then test three different cases for a\ = 25% and again three different cases for a\ = 55%. In 
all of these cases the condition of hitting the barrier comes into play and we can see the advantage 
of conditional sampling clearly. The QMC+LT method performs worse than MC+CS on one occasion 
but delivers good results for all other cases. Our new QMC+LT+CS method always outperforms the 
QMC+LT method and outperforms the MC+CS method by factors ranging from 1.6 to 8 with respect 
to the standard deviation. 

In Figure 1 we show the convergence graphs for the second and fifth examples. The convergence 
factors a, which we obtained from the linear regression log(c) = j3 — a\og(N), are shown in Table 2. 
We see that in two out of four cases our conditional scheme improves the convergence, and in the other 
cases the convergence is about the same as the original QMC+LT method. Overall, the results of our 
conditional sampling scheme QMC+LT+CS are very satisfactory. 
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Figure 1: Convergence graphs corresponding to the single barrier Asian basket. 





QMC+LT+CS 


QMC+LT 


MC+CS 


(P 1; 0.25, 125, 70) 


0.6991 


0.5363 


0.4939 


(P 2 , 0.25, 125, 70) 


0.7296 


0.5321 


0.4932 


(Pi, 0.55, 105, 70) 


0.5526 


0.5838 


0.5046 


(P 2 , 0.55, 105, 70) 


0.6145 


0.5828 


0.5043 



Table 2: The convergence factors a obtained from the linear regression log(er) — f3 — alog(iV) for single 
barrier Asian basket. 



5.2 Double barrier binary Asian 

Consider the following option payoff on a single asset: 

g = l\ — y^S(ti) > 100 1 1 1 min S(U)>B L \ll. max SUA < B u \ . 

I m j [_i=l,...,m J [_i=l,...,m J 

The model parameters are 5(0) = 100, a = 30%, T — 3 months and r = 0%. The results are presented 
in Table 3. 

In the first three examples, the number of observation dates m is fixed at 60, and the barriers are 
chosen increasingly closer to S(0) . As the barriers are tightened the performance of both the QMC+LT 
and QMC+LT+CS method drop, but the QMC+LT+CS scheme clearly keeps the upper hand. In order 
to push the QMC+LT+CS method to the limit, we consider the extreme example of very tight barriers 
(B L , B u ) = (98, 102) and choose m = 2,3, 4. These results are shown in the last three rows of Table 3. 
Although our method gets a serious variance reduction compared to the original QMC+LT method, we 
also notice that for m = 4 our QMC+LT+CS scheme does not outperform the MC+CS method. To 
understand what is happening, we can plot the projection of the payoff on the first two dimensions u\ 
and U2- This is shown in Figure 2. We see that the LT algorithm is not able to separate the positive 
payoffs from the zero payoffs for increasing m (where, e.g., in [18] it is observed that having a good 
separation will make the QMC method more effective). A possible solution here might be to use a non- 
linear method. Another observation is that in case of a tight double barrier, it is not always possible to 
find a valid bound on u\. This is illustrated in Figure 3 where the path is drawn for different values of U\ 
for m = 4 for a specific choice of ui, ■ ■ ■ ,1*4. The red lines indicate the barriers at 98 and 102. Clearly, 
there is no possible choice for u\ such that the asset path stays between these barriers. While this is a 
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/ dL r>U\ 
\ t> , _t> J 


m 






yi 7 1UUU j 


fin 


IDoO/O 


IDoO /o 


(50,150) 


60 


1692% 


1180% 


(90,110) 


60 


292% 


118% 


(98, 102) 


2 


296% 


26% 


(98, 102) 


3 


100% 


13% 


(98, 102) 


4 


40% 


4% 



Table 3: Double barrier binary Asian. The reported numbers are the standard deviations of the MC+CS 
method divided by those of the QMC+LT and QMC+LT+CS methods. The MC+CS method uses 
163840 samples, while the QMC+LT and QMC+LT+CS methods use 4096 samples and 40 independent 
shifts. 



drawback of our method, it must be noted that a product with such barriers is not encountered in the 
market. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



(a) m = 2 (b) m = 3 (c) m = 4 

Figure 2: The projection of the double barrier binary Asian payoff on the first two dimensions u\ and 
u 2 . Blue circles indicate a zero payoff, red circles a non-zero payoff. 



5.3 Single barrier Asian 

We now construct an example for which the first m elements of the first column of A do not have the 
same sign and we thus need both T u and Yd for the QMC+LT+CS algorithm. In this example we take 
the following payoff on two assets: 

\ (Si(h) + S 2 (h) + S 1 (t 2 ) + S 2 (t 2 )) > lj I{min (Si(h), S^h)) < B} . 

The model parameters are a x = 40%, a 2 = 60%, r = 8%, Si(0) = 5 2 (0) = 1, h = |, t 2 = 1 and B = 1.1. 
Straightforward calculations show that an < and a 2 \ > if 



-<Jl —2(7! 



a 2 e' 



i/ 2 ' ^(e^-^/ 2 )/ 2 + e r - CT i/ 2 ) 



-73.68%, -71 



We take p — —72%. In this case, it is possible that > T u and our sample is wasted. This happened 
in our simulation about 50% of the time. Still we get a very nice result as the standard deviation of the 
MC+CS method divided by those of the QMC+LT and QMC+LT+CS methods based on 4096 samples 
and 40 shifts result in 384% and 399% respectively. A convergence plot is given in figure Figure 4(a) 
while in Figure 4(b) the effect of u\ on a sample path for which > T u is given. 
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Figure 3: The effect of U\ on a sample path for the double barrier binary Asian with very narrow barriers. 




10 2 10 3 10* 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

N t 

(a) Convergence graph. (b) Effect of ui on a sample path. 



Figure 4: Single barrier Asian example with mixed signs for the conditional sampling. 
5.4 Binary barrier 

As a final experiment consider the following option payoff on one asset: 

g = I |max S'(ti) < 105 j . 

The model parameters are 5(0) = 100, a — 30%, T = 3 months and r = 0%. We are essentially 
integrating a closed volume in the m-dimensional hypercube. The resulting convergence factors are 
given in Table 4. These factors are calculated using the linear regression log(cr) = f3 — alog(jV). As 
we can see, increasing m decreases the convergence speed for both QMC+LT and QMC+LT+CS. This 
example illustrates that conditional sampling results in a variance reduction, but not necessarily in an 
increase in convergence speed. One possible explanation is as given in Section 5.2, i.e., due to the 
barrier condition the LT algorithm is not capable any more of aligning the discontinuities with the axes. 
Furthermore, as illustrated in [18], it might be necessary to sample the boundary of the closed volume 
adaptively to increase the speed of convergence. 

6 Extension: Knock-in options 

In the previous section we have constructed the conditional sampling scheme for the LT algorithm for 
knock-out options, i.e., the value is set to zero when crossing the barrier(s). In this section, we will 
extend our algorithm to knock-in options where the option is worthless unless the barrier is crossed. 
Such an extension was also considered in [9] where the authors note "Dealing with knock-in options is 
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m 


QMC+LT 


QMC+LT+CS 


2 


0.7804 


1.1830 


3 


0.6820 


1.0383 


4 


0.6815 


0.7925 


5 


0.5894 


0.7794 


6 


0.5328 


0.6015 


60 


0.5454 


0.5040 



Table 4: The convergence factors a, which we obtained from the linear regression log(er) = j3 — a log (AT) 
for the binary barrier. 



not simple, but is possible if there is a known expression gi{Si) for the present value of a barrierless 
option, received at time U when the state vector is Si, whose payoff will be g(S m ) at time t m ." They 
then continue with "This is the case for sufficiently simple knock-in options" and give some examples. 
We will first discuss the method described in [9] and then propose an easier and more flexible scheme 
which can also be used when this condition is not fulfilled. Note that in [9] paths are constructed in an 
incremental fashion as this is key to understanding their procedure. 
We consider an up-&-in condition, i.e., a payoff 

g(S 1 (ti),...,S 1 (tm),---,S n (t m ))=max(f(S 1 (t 1 ),...,S n (t m )),0)l\ max S^tj) > B 

|_j=l,...,m 

We first explain the scheme from [9], Assume that there is a known expression for the option value 
without the barrier at every time step given the state vector Si = (S'i(ti), . . . , S%(ti), . . . , S n (ti)) (an 
incomplete path from to up to tj), denoted by g(Si), then by noting that a knocked-in barrier option 
becomes a regular option on knock-in, one samples the variable 

f;(l{. = nua_ i 5 1 (t i )<B}-l{.max i 5 1 (t j )<B]] 5 (5 i )= £ g{Si). 

i=l ^ ^ * ^ i=argmin 3 = 1 m Si(tj)>S 

To take advantage of importance sampling the authors in [9] extend this basic sampling scheme so that 
in each time step tj, there are two successors to Si(fj): Si(U + i) simulated conditional on no knock- in 
and £*(tj+i) simulated conditional on knock- in. (The incremental path is afterwards continued from the 
Si{ti+i) value.) The authors then propose to sample 

m 

Li-i (1 - P [Si(ti-i) < B | Sx(tj)}) g(S*), 

t=i 

where Li-\ is the likelihood as defined in (4). Now each term in the sum contributes to the random 
variable. However, this scheme is rather restrictive as it is required to know the value of the option in 
each time step. When considering a knock-in arithmetic Asian option for example, this is not the case. 

We propose an easier and more flexible scheme as an extension of our method from Section 4. Under 
the LT algorithm, when considering an up-&-in option, we must have for some tj 

o-i Wiitj) > b{t 3 ), 

where b(tj) is defined in (6). The difference with a knock-out type option is that now the above condition 
must only hold for at least one tj. Again suppose aj i > for j e V C {1, . . . ,m} and dj t i < for all 
other j ^ V . For an up-&-in option we can essentially "force" the paths to cross the barrier level by 
imposing 

b(tj) - a jt 2 i S>~ 1 (u2) - ... - a ] , mn '&~ 1 {um,n) 

a i,i 
or 

b(tj) - aj > 2^ > ~ 1 (^2) - ■■■ - aj, mn ^ 1 ( u mn) ' 



$ (tti) > min 



$ 1 (ui) < max 



14 




Figure 5: Verification with continuous down-&-in put option with knock-in condition. 



For a down-&-in option the condition on u\ is analogously 

b(tj) - aj^fl? -1 ^) ^ ■■■ - a^mn^^iUmn) 

or 

b(tj) - aj,2^ > ~ 1 (M 2 ) - ... - aj-mn^Vmn) " 

We will denote these bounds again by and T u . We can again prove that this estimator is unbiased 
and has a standard deviation at most that of the unaltered QMC+LT method. 

Theorem 6. The estimator for knock-in options based on sampling 

92 = (T„ - T d ) max(/(ui, u 2 , u mn ),Q) (8) 
is unbiased. That is, E[g] = £[52]- 

Proof. Analogous to Theorem 3. □ 

Theorem 7. When using a randomly shifted quasi-Monte Carlo method, the estimator defined in (8) 
has reduced variance. That is, Var[g 2 ] < Vax[g]. The inequality is strict if ¥ [maxj Si(tj) > B] < 1 and 
E [/] > 0, i.e., if there is a positive chance of no knock-in and a positive payoff. 

Proof. Analogous to Theorem 5. □ 



$ 1 (ui) < max 



$ 1 (ui) > min 



6.1 Verification by continuous down-&-in put with knock-in condition 

We can check the valuation of a simple continuous down-&-in put option to illustrate the unbiasedness 
of our method for pricing knock-in options. Consider the put variant: 

9(S(t)te[o,T]) = max (K - S{T),0) I / min S(t) < B j . 

We approximate the option value by using a fine time discretization of 500 steps. The value of the option 
using a continuous barrier is known analytically, see, e.g., [24]. The model parameters are K = 100, 
B = 80, g = 20% and r — 5%. Three maturities are considered: T = 1 year, T = 6 months and T = 1 
month. The valuation is based on 4096 samples and 40 independent digital shifts. The result is shown 
in Figure 5. We see that the valuations are very close to the analytic values. 
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(P rr-, P. K ) 




(P rr, R K\ 




( Pi n 25 1 05 70") 


287% 


( Pn 25 1 05 701 


237% 


(Pi, 0.25, 125, 70) 


556% 


(P 2 ,0.25, 125,70) 


320% 


(Pi, 0.25, 200, 70) 


2272% 


(P 2 , 0.25, 200, 70) 


345% 


(Pi, 0.55, 125, 70) 


297% 


(P 2 , 0.55, 125, 70) 


383% 


(Pi, 0.55, 125, 90) 


213% 


(P 2 , 0.55, 125, 90) 


372% 


(Pi, 0.25, 120, 110) 


180% 


(P 2 , 0.25, 120, 110) 


178% 


(Pi, 0.25, 110, 100) 


133% 


(P 2 , 0.25, 110, 100) 


171% 



Table 5: Basket Asian with knock-in condition. The reported numbers are the standard deviations of 
the QMC+LT method (as the MC+CS method does not work for this example) divided by those of the 
QMC+LT+CS method. Both methods use 4096 samples and 40 independent shifts. 



6.2 Basket Asian with knock-in condition 

As a second test we consider the same Asian basket on four assets as in Section 5.1 but now with a 
knock-in condition. Since there does not exist a closed-form solution of the arithmetic Asian option, 
we are not able to use the MC+CS method from [9]. Therefore, we report the ratios of the standard 
deviations of the value estimates of QMC+LT to QMC+LT+CS for various choices of oi, K and B in 
Table 5. The first three examples take a x = 25%, K = 70 and B = 105, 125 and 200. When the barrier 
is close to the starting value of Si (being 100), the probability of knocking- in is higher, but we still see a 
significant variance reduction when using conditional sampling. As the barrier moves further away from 
the initial asset value, the probability of hitting the barrier drops, and we see that the variance reduction 
improves even further. 

The next two examples take U\ = 55%, B = 125 and K = 70 or K — 90. When K = 70, we see 
a variance reduction for the QMC+LT+CS method, but we see that increasing <7i from 25% to 55% 
decreases the variance reduction for Pi. This is what we would expect, since a higher volatility for Si 
implies that more paths are hitting the barrier without conditional sampling. 

Also for the last two examples we observe that QMC+LT+CS shows a variance reduction to the 
QMC+LT method. In conclusion: the results of our conditional sampling scheme QMC+LT+CS exten- 
ded for knock-in options are again very satisfactory. 



7 Extension: Root-finding 

In this section we want to exploit the fact that our construction in Section 4 only modifies the first 
uniform variate to satisfy the barrier condition. Here we propose to exploit the influence this variate 
has by calculating the bounds on ui as in Section 4 to fulfill the barrier condition and then, using root- 
finding, determine the region (between these bounds) which produces a positive payoff. We can then 
integrate out ui analytically. 

The idea of using root-finding is quite natural and has been previously discussed in, e.g., [6] and 
[10], albeit more directly. For this method to be applicable it is necessary that we are able to determine 
bounds on ui (or equivalently zi = $(iti)) for it to be applicable. We will consider an Asian basket 
as an example of a more complicated payoff and then afterwards a simple put option with a knock-in 
condition. 



7.1 Single barrier Asian basket 

An Asian basket option on n assets has a payoff 

C. n m \ 
— E S ll (U 2 ) ~ K,0 !{•}, 
ran * — ' * — ' / 

with possible some barrier condition for knock-out or knock-in (denoted by !{•} above). The following 
theorem ensures that the problem is well-posed for Asian barrier options. 
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Theorem 8. Given any covariance matrix £ with factorization A A 1 , the function 

1 mn 

fizA = — y > 5 , il (o) e ( r - CT< i' /2 ^2+5:r=?2"«^ e Q * 121 -k 

mn ^— ' 

i=l 

{where i\ = \_(i—l)/m\ +1 arid 12 = i — (ii — l)m) /ias most two zeros. If all elements an, i = 1, . . . , mn 
are of the same sign and K > 0, t/ien the function f{z\) has exactly one zero. 

Proof. First, consider the case where an > for all i. Then each asset price is a monotonically increasing 
function of z\. Therefore, the sum is a monotone function of z\ as well. Since Si t (ti 2 ) — > for z\ — > — oo 
and <S , i 1 (tj 2 ) — > +oo for z\ — > +oo, there is exactly only one zero of the function f{z\) if K > 0, and 
none if K < 0. A similar argument can be made when an < for all i. 
Consider then the case where an has mixed signs. We have that 

dzi mn * — ' 

1 i=l 

and 

,2 h mn 

f'(zi) = = — y / S n (0)e {r -^ /2)t ^ + ^ a ' kZk a 2 tl e a ^ . (10) 

1 i—l 

Clearly, f"(zi) > everywhere. The function f'{z\) can be written as the difference of a monotone 
increasing function and a monotone decreasing function in Z\. This implies that f[z\) — > — oo as 
z\ — > — oo and f'{zi) — > +oo as z\ — > +oo. So f'{zi) changes sign at least once. For f(z\) to have more 
than two zeros, /'(zi) would have to change sign more than twice. However, f'{zi) has only one zero, 
which can be easily seen using Rolle's theorem and the observation that f"{z\) > 0. Therefore f'(zi) 
changes sign exactly once, which leaves the possibility of zero, one or two zeros. □ 

In case the elements in the first column of A have the same sign, we are guaranteed there is exactly 
one zero of the payoff function, which can be found using a root-finding algorithm such as Newton- 
Raphson. In case the root falls outside the interval [T^, T n ] as defined in Section 4, there is no value for 
z\ which produces a positive payoff and we set the sample value equal to zero. If, on the other hand, 
the root does fall inside the interval, we can analytically integrate the payoff function over z\ (or u\) by 
using the forthcoming lemma. 

When the elements in the first column of A have mixed signs, things get a bit more complicated since, 
given Z2 to z mni there might be zero, one or two roots of the payoff function g. We do know however 
that /' has exactly one zero, so we propose to find this zero using an appropriate root-finding algorithm, 
and then to check whether / is positive, zero or negative in this point. In case / is positive, there is no 
zero and we can use the bounds found for the barrier condition as the range of integration. We can do 
the same if / is zero. If / is negative, we must find two zeros, again using the appropriate root-finding 
algorithm, initiated with a point to the left and to the right of the zero of /'. 

If we call these roots Td and T u , where they can be respectively equal to — oo and oo, then we can 
define = max(T ( j,r ( j) and S„ = min(Y n ,r„) an integrate out the first dimension analytically as the 
next lemma shows. 

Lemma 1. Given Z2 to z mn and the bounds Srf(z2, . . . , z mn ) and S„(z2, . . . , z mn ) on z\ for which an 
Asian barrier option has a positive payoff and satisfies the barrier condition, then the mean value over 
Zi is given by 

E*i[ff(*i)l ^2 1 ■ ■ ■ i Zmn\ 

1 mn 

= — V S ll {Q)e^-^' 2 ^+^ «<^ e <4/2 ($( Su _ aa ) _ _ aa )) - ($(£„) - $(Z d ))K. 

mn. * — * 
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^lvl\_^-t-ljl -fuO 


(Pi 25 125 10(11 


2039% 


958% 


c Po n 25 125 inn) 


1 955% 


1172% 


(Pi, 0.25, 110, 100) 


960% 


446% 


(P 2 , 0.25, 110,100) 


1172% 


595% 


(Pi, 0.25, 105,100) 


638% 


263% 


(P 2 ,0.25, 105,100) 


761% 


367% 


(Pi, 0.25, 110,90) 


910% 


737% 


(P 2 , 0.25, 110, 90) 


908% 


782% 


(Pi, 0.25, 105,90) 


757% 


576% 


(P 2 , 0.25, 105, 90) 


623% 


581% 


(Pi, 0.25, 125,110) 


1923% 


489% 


(P 2 , 0.25, 125,110) 


4693% 


683% 


(Pi, 0.55, 125,100) 


939% 


437% 


(P 2 ,0.55, 125,100) 


1082% 




(Pi, 0.55, 125, 110) 


1035% 


234% 


(P 2 , 0.55, 125,110) 


1381% 





Table 6: Single barrier Asian basket. The reported numbers are the standard deviations of the MC+CS 
method divided by those of the QMC+LT+CS and QMC+LT+CS+RF methods. The MC+CS method 
uses 163840 samples, while the QMC+LT+CS and QMC+LT+CS+RF methods use 4096 samples and 
40 independent shifts. 



Proof. We can write 



mn 



E zl [g{ Zl )\z2, . . . , z mn ] = V 5^(0) e < r - CT n/^+STIW** 



i=l 



f3 u e anx-x /2 
5, 



dx-K 



-x 2 /2 



-dx. 



The result then follows by straightforward calculations. □ 

The extra cost of using the root-finding procedure is limited: in the worst case (i.e., when the elements 
in the first column of A have mixed signs) we need to calculate (9) and (10). We can actually precomputc 

which has to be calculated in the original algorithm as well, and then use vector multiplication with 
(al 1 e ail )i, where the power of an corresponds to the rth derivative of /. The original algorithm re- 
quired one such multiplication, so if we take for example four steps for our root-finder, we need three 
multiplications on top of the original algorithm. Higher derivatives can be computed at almost no ex- 
tra cost, since the only thing that changes is the power r. In our algorithm, we use a fourth-order 
Newton-Raphson algorithm with 10 iteration steps. 

At first thought one might think that having a negative correlation might imply different signs in the 
first column of A. Surprisingly, as shown in Section 5.3 this relation does not hold necessarily. Three 
examples of the root-finding problem are shown in Figure 6, where we consider an Asian up-&-out option 
with parameters S"i(0) = 5 2 (0) = 100, r = 5%, a 1 = 30%, a 2 = 40%, B = 110, K = 100, T = 1, m = 260 
and different correlations p between the two assets to display some possible cases for the roots of f{z\). 
Figure 6(b) is an example with negative correlation p — —30% but only one zero. 




Figure 6: Possible root configurations for an Asian basket on two stocks with different correlation. 

For our numerical test we look again at the single barrier Asian basket from Section 5.1. The results 
for different choices of K, B and o\ are shown in Table 6. By QMC+LT+CS+RF we denote the new 
algorithm from this section. There is a significant variance reduction compared to the other methods for 
all parameter choices. 
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(a, B, K) 


QMC+LT+CS+RF 


QMC+LT+CS 


(0.25,90,100) 


442% 


184% 


(0.25,90,90) 


551% 


134% 


(0.25,90,80) 


742% 


176% 


(0.25,75,80) 


884% 


480% 


(0.55,75,80) 


575% 


169% 


(0.55,70,60) 


657% 


185% 



Table 7: Down-&-in put option. The reported numbers are the standard deviations of the QMC+LT 
method divided by those of the QMC+LT+CS and QMC+LT+CS+RF methods. All methods use 4096 
samples and 40 independent shifts. 



7.2 Down-&-in Put option 

As a final example we consider a simple put option with a knock-in condition: 

g = max (if - S(t m ), 0)1 { min S(U) < b\ . 

We fix m = 130, S(0) = 100, r = 5% and T — 6m. Table 7 shows the results for different values of a, B 
and K. In this example of a knock-in option, we are able to force 100% of the samples to result in positive 
payoffs. For a wide variety of parameter choices, including different volatilities and both in-the-money 
and out-the-money barriers and options we see significant gains compared to the QMC+LT+CS and 
QMC+LT methods. 

8 Conclusion 

In this paper we have devised a conditional sampling method for QMC sampling under the LT algorithm 
for barrier options. Our method tries to satisfy the barrier condition by modifying the first uniform 
variable used to construct the sample paths. We have shown that the method is unbiased for randomized 
QMC methods and has a variance which is lower than or equal to the method without using conditional 
sampling. Furthermore, the method can be used for both knock-out and knock-in conditions with the 
same basic principle of modifying the first uniform variable. When adding a root-finding method we can 
additionally also satisfy to have always a positive payoff which further improves the performance of the 
method. Extensive numerical results show the effectiveness of the new method. 

References 

[1] R. Cools, F. Y. Kuo, and D. Nuyens. Constructing embedded lattice rules for multivariate integra- 
tion. SIAM Journal on Scientific Computing, 28(6):2162-2188, 2006. 

[2] E. Derman and I. Kani. The ins and outs of barrier options: Part 1. Derivatives Quarterly, 2(3):55- 
67, 1996. 

[3] E. Derman and I. Kani. The ins and outs of barrier options: Part 2. Derivatives Quarterly, 3(3) :73- 
80, 1997. 

[4] J. Dick and F. Pillichshammer. Digital Nets and Sequences: Discrepancy Theory and Quasi-Monte 
Carlo Integration. Cambridge University Press, 2010. 

[5] C. Fries. Mathematical Finance: Theory, Modeling, Implementation. Wiley, 2007. 

[6] T. Gerstner. Sparse grid quadrature methods for computational finance. Habilitation thesis, 2007. 

[7] M. B. Giles, F. Y. Kuo, I. H. Sloan, and B. J. Waterhouse. Quasi-Monte Carlo for finance applica- 
tions. ANZIAM Journal, 50:308-323, 2008. 



19 



[8] P. Glasserman. Monte Carlo Methods in Financial Engineering. Springer, 2003. 

[9] P. Glasserman and J. Staum. Conditioning on one-step survival for barrier option simulations. 
Operations Research, 49(6):923-937, 2001. 

[10] M. Holtz. Sparse Grid Quadrature in High Dimensions with Applications in Finance and Insurance. 
Springer- Verlag, 2011. 

[11] J. C. Hull. Options, Futures and Other Derivatives (6th Edition). Prentice Hall, 2005. 

[12] J. Imai and K. S. Tan. A general dimension reduction technique for derivative pricing. Journal of 
Computational Finance, 10(2):129-155, 2006. 

[13] S. Joe and F. Y. Kuo. Constructing Sobol' sequences with better two-dimensional projections. SIAM 
Journal of Scientific Computing, 30:2635-2654, 2008. 

[14] M. Joshi and R. Tang. Pricing and deltas of discretely-monitored barrier options using stratified 
sampling on the hitting-times to the barrier. International Journal of Theoretical and Applied 
Finance (IJTAF), 13(5):717-750, 2010. 

[15] C. Lemieux. Monte Carlo and Quasi-Monte Carlo Sampling. Springer Series in Statistics. Springer, 
2009. 

[16] C. Lemieux and P. L'Ecuyer. Variance reduction via lattice rules. Management Science, 46(9):1214- 
1235, 2000. 

[17] C. Lemieux and P. L'Ecuyer. Randomized polynomial lattice rules for multivariate integration and 
simulation. SIAM Journal on Scientific Computing, 24(5): 1768-1789, 2002. 

[18] D. Nuyens and B. J. Waterhouse. A global adaptive quasi-Monte Carlo algorithm for functions of 
low truncation dimension applied to problems from finance. In H. Wozniakowski and L. Plaskota, 
editors, Monte Carlo and Quasi-Monte Carlo Methods 2010, pages 1-18. Springer- Verlag, 2012. To 
appear. 

[19] A. Papageorgiou. The Brownian bridge does not offer a consistent advantage in quasi-Monte Carlo 
integration. Journal of Complexity, 18(1):171-186, 2002. 

[20] I. H. Sloan and S. Joe. Lattice Methods for Multiple Integration. Oxford Science Publications, 1994. 

[21] J. Traub and A. Papageorgiou. Beating Monte Carlo. Risk, pages 63-65, 1996. 

[22] J. Traub and S. Paskov. Faster valuation of financial derivatives. Journal of Portfolio Management, 
22(1):113-120, 1995. 

[23] X. Wang and I. H. Sloan. Quasi-Monte Carlo methods in financial engineering: An equivalence 
principle and dimension reduction. Operations Research, 59(l):80-95, 2011. 

[24] P. Wilmott. Paul Wilmott on Quantitative Finance. Wiley, 2006. 

[25] U. Wystup. FX options and structured products. Wiley, 2006. 

[26] R. Zvan, P. Forsyth, and K. Vetzal. Discrete Asian barrier options. Journal of Computational 
Finance, 3:41-68, 1999. 



20 



